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We present a variational method to calculate the exchange interactions among donor clusters in a 
semiconductor. Such clusters are candidates for a so-called control-qubit architecture for quantum 
information, where the effective exchange coupling between two atoms is controlled by the electronic 
state of a third. We use a combination of the effective-mass approximation and the quantum defect 
method; our variational ansatz is particularly suited to cases where an excited state of one of the 
donors (control) is partially delocalised over several different centres, forming an analogue of an 
extended molecular orbital. Our method allows calculations of the "on/ off' ratios of exchange 
interactions in such cases. We compare exchange interactions when the control is in the "on" 
and "off' states, and find that both the magnitude and sign of the exchange interactions may be 
changed. To rationalize the sign-change, we carry out a simple Green's function perturbation-theory 
calculation. This simple model qualitatively explains the sign change and illustrates its origins both 
in ring-exchange processes and in the delocalization of the control electron. We also compute 
probability distributions for the coupling strengths over the ensemble of clusters, and show that 
excitation of the control causes narrowing of the distributions along with shifts to larger magnitudes 
and from anti-ferromagnetic to ferromagnetic coupling. 

PACS numbers: 03.67.Lx, 71.55.Cn 
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I. INTRODUCTION 

Proposed silicon-based quantum computer architec- 
tures [TJ [5] have attracted attention because of their 
promise of scalability and their potential for integration 
with existing CMOS architectures. Localized spins in Si 
are one very promising means to represent quantum in- 
formation: spin-i objects are natural two- level systems, 
so they may naturally be used to embody a qubit. Fur- 
thermore, in pure Si, electron spins are associated only 
with defects and so are naturally isolated. 

Implementing general quantum gates requires one to 
produce entangling interactions between localized spins 
in order to evolve these spin states. Previous ap- 
proaches to controlling these interactions [T] required 
gate electrodes positioned near to specific highly polar- 
izable (therefore shallow) donors; such defects are how- 
ever readily ionized except at low temperatures, and have 
spin-lattice relaxation times declining rapidly with tem- 
perature because of the presence of spin- flip Raman tran- 
sitions P] . In addition the presence of the gates may in- 
troduce significant additional decoherence through their 
interaction with the polarizable defects, as well as posing 
formidable difficulties in the accurate positioning of the 
defects relative to the electrodes [H |5] . 

An alternative scheme for controlling the interactions, 
which avoids the need for electrodes, was proposed in [6]. 
Quantum bits are encoded into electron spins of deep 
donors, with the advantages of lower ionization prob- 
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abiliities and much longer spin-lattice relaxation times 
[S]. Interactions are controlled, not by gates, but by lo- 
cal electronic excitations f.6|. Specifically, let (A,B) be 
two such deep donor atoms. Their spacing should be suf- 
ficiently large that the ground-state interaction between 
donor spins is small (ideally negligible). Controlled op- 
tical excitation [7l |8] can promote a "control" electron 
from a nearby impurity C into an excited state that is to 
some extent delocalized across A, B and C. In this excited 
state, there is expected to be an effective exchange inter- 
action induced between the two qubit spins. Qubit-qubit 
interactions are therefore switched on by optical excita- 
tion and off by (stimulated) de-excitation of the control 
electron. Note that in this picture the electronic excita- 
tion is real, not virtual, in contrast to other schemes that 
have been proposed to couple quantum dots in semicon- 
ductors 7; this has the advantage that the wavelength 
variations arising from the inhomogeneity of the sample 
can be exploited to address it on scales much smaller than 
the optical wavelength [B^ , while leaving the control spins 
unentangled with the qubits if the operating parameters 
are chosen carefully [10]. 

To build a quantitative model of these processes, we 
need to understand how the exchange interactions de- 
pend on the different types of donors involved (deep or 
shallow), their separations, the presence or otherwise of 
electronic excitations, and the extent of delocalization of 
these excitations. Exchange interactions between pairs of 
hydrogenic impurities in semiconductors have been stud- 
ied for many years [HJ [H [H [H [TU [H] , and we re- 
cently showed how these approaches could be easily gen- 
eralized to pairs of deep donors with an arbitrary bind- 
ing energy |17j using a simple quantum-defect approach. 



The exchange among defect clusters has been much less 
studied, except where the clusters strongly influence the 
bulk magnetic properties, as in for example (Ga,Mn)As 
[18j , although the general importance of multi-centre in- 
teractions in determining exchange properties has been 
known for many years, notably in the context of solid '^He 
[19| . Their potential significance for quantum informa- 
tion processing with electron spins was pointed out more 
recently ^U\ [3T1 [52] ; however this work focussed on har- 
monic potential wells in highly symmetric arrangements. 
Deviations from the classic form of Heisenberg exchange 
have also been discussed [23] in the context of multiple 
quantum dots, and are important when the dots are very 
strongly coupled or in very large magnetic fields. 

We focus on a different aspect of the problem, and 
present a simple variational calculation designed to cap- 
ture the essential physics: the presence of deep and shal- 
low impurities, a range of geometries, and the possibility 
of electron delocalization in the excited state. We adopt 
a reference model in which the qubits are deep donors 
(for example, Bi:Si) while the control atom is a shallow 
donor (e.g. P:Si), and treat it in configurations where the 
control is equidistant from the qubits. This remaining 
symmetry simplifies the calculations while enabling us to 
study separately the effects of qubit-qubit and control- 
qubit separations. Our basic model can be thought of 
as a multi-spin generalization of the well-known Heitler- 
London approach [531 123 to the exchange interactions 
in the hydrogen molecule. Although this approximation 
fails at very large separations [^ because of the neglect 
of correlation terms, it is very accurate over the phys- 
ically interesting range of separations up to about 12 
effective Bohr radii, beyond which dipolar interactions 
anyway become dominant |17j . We also show that we 
can rationalize our main results using simple perturba- 
tion theory. 

The remainder of the paper is organized as follows. In 
j plj we introduce our methods and the basis used to de- 
scribe the wave functions of deep and shallow donors. In 
the §III| we present our results for a range of qubit and 
control positions. In pV| we discuss the results and show 
how they can be related to the so-called kinetic exchange 
arising from multi-centre interactions by using perturba- 
tion theory; then in Sjv] we show how the distributions 
of exchange that would be probed by bulk measurements 
are influenced by the properties of the defect clusters and 
their degree of excitation. We do not address the opti- 
cal excitation and de-excitation processes; these will be 
discussed in a forthcoming paper |26j . 



II. METHOD OF CALCULATION 
A. Geometry and representation of the solid 

Let Re, Ra, and Rb be the position vectors of the con- 
trol nucleus, the flrst qubit nucleus, and the second qubit 
nucleus respectively. We take the distance between the 
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FIG. 1: The geometry used for the three-atom system. 



qubits as \Ra — Rb\ — Dqq] for reasons that will become 
clear below, we take the triangle formed by our three 
atoms to be isosceles, so that a reflection plane exists bi- 
secting the two qubit atoms. Let the height of this trian- 
gle be H] the control atom is then equidistant from both 

qubits and \Rc - Ra\ ^ \Rc - Rb\ = y/lP+jD/W = 
DcQ ■ This geometry is sketched in Figure [Tj 

We treat the three donors within isotropic, single-band 
effective-mass theory [571 HH] , introducing quantum de- 
fect corrections [HU] where necessary to modify the bind- 
ing energies. This is the simplest model that enables 
us to concentrate on the mutual interaction of the three 
donors, and we have recently shown [T7] that it is able to 
describe the exchange between donor pairs with arbitrary 
binding energy. However, it does neglect the valley de- 
generacy of the silicon conduction band, which in reality 
contributes a degeneracy of the donor states and corre- 
sponding interference terms in the exchange [12l |13l [14] . 
In common with other recent treatments we do not treat 
the conduction-band anisotropy directly, instead captur- 
ing its main effects by choosing an appropriate average 
effective mass. 

Within the effective-mass model, the Hamiltonian of 
this three-electron system is 



H 



E^^nE 



1 



where 



Zc 



2 ^. \r. 



Za 



Zb 



Re 
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(1) 



(2) 



Here r^ labels the coordinate of electron i, and Zc, Za, 
and Zb are respectively the charges of the control nu- 
cleus, the first nucleus, and the second nucleus. We con- 
sider single donors and treat only the outermost electron 
explicitly and therefore take Zc = Za = Zb = 1. Since 
we neglect spin-orbit coupling and dipole-dipole interac- 
tions, the system Hamiltonian contains no explicit spin 
operators and the total spin S' is a good quantum num- 
ber. 

We use scaled atomic units throughout: we take a 
mean (directionally averaged) conduction-band effective 



mass for Si of m* = O.SSrrie (where rrie is the free elec- 
tron mass), and a relative permittivity Sr = 11.4, leading 
to an effective Bohr radius about Og = 1.94 nm and an 
effective Hartree Ha* = 33.48 meV JW\- 



B. Single-centre basis functions 

We suppose there is a single important orbital (j)c at 
the control site (whose nature changes depending on the 
degree of excitation of the control) and two important or- 
bitals at each qubit site: ground-state orbitals Xi (where 
i G {A, B}) which accommodate the qubit electrons, and 
excited orbitals (j>i which allow for delocalization of the 
control electron. In this paper we will assume that all 
of these states have s symmetry, and use quantum-defect 
methods to represent their radial parts. 

We make this choice so that we can work equally well 
with deep and shallow impurities (deep donors are espe- 
cially promising for QIP applications [B]), and still have 
functions that are described by effective mass theory far 
from the defects. Bebb ^U\ showed how quantum defect 
methods could be used to obtain good approximate wave 
functions in the region outside the impurity ion core, us- 
ing only a knowledge of the energy eigenvalues. In this 
region, where the potential approaches u{r) ~ — 1/r, the 
radial Schrodinger equation is 



2dr2 



2r2 



1 



-uir)-^)P{r) = 0, 



(3) 



where — l/2i/2 is just e, the observed energy level. It is 
conventional to write i^ = n — fi(en), where n is the usual 
principal quantum number, fi is the quantum defect, and 



Ha* 
2^ 



Ha* 



2(n-Ai(e„))2 



(4) 



where Ha* is the effective Hartree. 

The solution of ^ that is regular at infinity is P^_i{r), 
but if /i 7^ this solution diverges at the origin if the 



potential remains 1/r down to short distances. We must 
therefore suppose that the short-range deviation of the 
potential from the Coulomb form is such that the true so- 
lution is indeed regular at the origin; the solution Pij,i{r) 
is only valid outside the core. In principle one could 
determine the entire eigenfunction using some form of 
central cell correction f3TJ |3H [M| • However for our pur- 
poses the asymptotic form is sufficient to calculate the ex- 
change interactions when the donor separations are large. 
The normalized quantum defect wavefunction is 



Pu,{r) = N,,iW,,..{2rlv). 



(5) 



Here W is a, Whittaker function ^B] , defined in terms of 
hypergeometric functions U by 



w^Az) 



g-./2^1/2+Mj;( 



1 



n + 11,1 + 2^1, z). (6) 



The approximate normalization constant N{v, I) is given 
by [30,^ 

In this paper we take ^ — Q (and hence v — 1,2,...) 
for shallow (control) donors, and /i = 0.3 (hence v = 
0.7, 1.7 . . .) for the more deeply bound qubit states. 



C. Variational method 

We make a simple variational assumption for the wave 
function of the control electron, writing the correspond- 
ing one-electron state as 



■0c = a0c + P(t>A + l<t>B- 



(7) 



Although this is a very simple choice, and will not reflect 
the full complexity of the problem, it does give the con- 
trol electron variational freedom to adjust between the 
extremes of being localized entirely on the control atom, 
or entirely on a "molecular" excited state of the qubit 
atoms A and B. 

We then construct a three-electron variational function 
in the form 



*totai(2^i'2;2,a;3) = A[i^c{'ri){XA{r2)XB{rz) - XA{r3)XB{r2))Q] 
A[^pc{ri){XA{r2)XB{r?,) - XA{rz)XB{r2))Di] 
A[^c{ri){XA{r2)XB{n) + XA{rz)XB{r2))DQ] 



*totai(a:i'a;2,a;3) 



*toL(2^i'a;2,a;3) 



(8) 

(9) 

(10) 



where Q refers to the quartet (total spin S = 3/2), Di 
refers the doublet (total spin S = 1/2) constructed from 
the triplet state of the qubits, and Dq refers to the dou- 
blet (total spin S — 1/2) constructed from the singlet 
state of the qubits. A is the antisymmetrizing operator. 
The existence of a plane of symmetry in our assumed ge- 



ometry means that the functions ^ total ''^^'^ ^ totals which 
have opposite parities with respect to the exchange of the 
qubits, cannot mix and can be varied independently. 

Note that, in constructing the above function, we have 
neglected high-energy 'ionic' configurations in which both 
qubit electrons reside on the same atom. We have also 
neglected correlations between the control electron and 



the qubit electrons. 

We now perform a variational calculation using unde- 
termined multiplier methods to minimize the expectation 
values of the system Hamiltonian: 



dE, ^^dS, 
da da 


==0 


(11) 


dE, , dS, 
dp^^dp 


= 


(12) 


dE, ^^dS, 
57 d'-f 


= 0, 


(13) 
(14) 



where i G {Q, Di, -Do}, -^ is the undetermined multiplier, 
K; = (^,;|iJ|*,;) and 5,; = (*,,|*,:). 



D. Extraction of the spin Hamiltonian 

Once we have the energies {Eq , Ejj-^ , Eu^ ) , of the spin 
states, we can calculate the exchange constants. Pro- 
vided spin-orbit coupling is negligible, the three-electron 
spin Hamiltonian in zero external field must take the form 

Heff == JCQASA ■ Sc + JcQBSb ' SC + JqQSA ' 83, (15) 

because this is the only form which is (i) invariant un- 
der simultaneous rotation of all the spins and (ii) even 
under time-reversal. Furthermore, atoms A and B are 
equivalent in our calculation, so Jcqa = Jcqb = JcQ- 
We can then calculate the exchange constants in the ef- 
fective Hamiltonian from the eigenvalues and eigenstates. 
As explained in piC| and Appendix]^ the doublet states 
are distinguished by the 'intermediate' angular momen- 
tum of the qubit spins (A and B); once these states have 
been identified, the parameters in equation ( 15 1 can be 
extracted to give 



Jqq - 

JcQ — 2\Ef 



Edo 



(16) 



where the zero of energy is chosen to be the weighted 
mean of the spin energies, so the spin Hamiltonian is 
traceless. Note that, with the sign convention adopted in 
setting up the spin Hamiltonian (151, positive couplings 



J correspond to antiferromagnetic interactions. 



III. VARIATIONAL CALCULATION: RESULTS 

In Figures [2] and |3] we show respectively the control- 
qubit coupling Jcq and the qubit-qubit coupling Jqq, 
as the qubit-qubit distance Dqq varies from 1 to 12 Og 
and the height H of the triangle formed by the qubits 
and control varies from to 10 Oq. In these and subse- 
quent figures the geomtry is as shown in Figure 111 the 
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FIG. 2: (Colour online.) The control-qubit exchange coupling 
Jcq (logarithmic scale) for the case where a shallow-donor 
control and deep-donor qubits {v — 0.7) are in the spatial 
configuration shown in Figure [l] as a function of qubit donor 
distance Dqq and the qubit-control distance Dcq- (a) Con- 
trol is in the ground state {u = 1.0); (b) control in the excited 
state {v — 2.0). The sign of exchange coupling is coded in 
color; red: positive (antiferromagnetic) and green: negative 
(ferromagnetic). Dqq varies from 1.0 to 12. Oa^; H varies 
from to lOaJ). Note that the vertical scales in the two plots 
are the same. 



two qubits are deep donors, and the shallow control is ci- 
ther (a) in the ground state {f = 1) or (b) in the excited 
state {ly — 2). All results are shown as a function of Dqq 

and the qubit-control distance Dcq — \ I H'^ 



dIqI^; 

the constraint Dqq > Dqq/ 2 is therefore enforced by 
our geometry. 

It is immediately evident from Figure 0^a) that even 
with the control in its ground state, Jcq depends 
(weakly) on Dqq as well as (strongly) on Dcq- Through- 
out most of the region Jcq is positive (antiferromag- 
netic), as would be expected for a two-electron system 
[34' , although in some small regions it becomes negative 
(ferromagnetic). Figure pTb), where the control is in the 
excited state, shows very different behaviour: now the 
effective coupling is much larger and decays much more 
slowly (as we would expect) but is also ferromagnetic 
over a significant part of parameter space (which is more 
surprising) . We will see below that this change is related 
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FIG. 3: (Colour online.) The qubit-qubit exchange coupling 
Jqq (logarithmic scale) for the case where a shallow-donor 
control and deep-donor qubits {v — 0.7) are in the spatial 
configuration shown in Figure [l] as a function of qubit donor 
distance Dqq and qubit-control distance Dcq- (a) control 
is in the ground state [v — 1.0); (b) control in the excited 
state {i' — 2.0). The sign of exchange coupling is coded in 
color; red: positive (antiferromagnetic) and green: negative 
(ferromagnetic). Dqq varies from 1.0 to 12.0oo; H varies 
from to IOoq . Note that the vertical scales in the two plots 
are the same, and that the constraint Dcq > Dqq is enforced 
by our geometry. 



to a change in the spatial location of the control electron, 
as it becomes delocalized over the three centres. 

Similar cooperative behaviour is evident in the be- 
haviour of the qubit-qubit coupling Jqq, which displays 
a more complex structure than the exponential decrease 
with Dqq that one might naively expect. Figure ^a) 
shows that the coupling is predominantly exponential 
in Dqq only when Dqq <^ Dqq, but this behaviour 
changes when the control atom becomes situated close 
to the midpoint of AB. For Dqq > Dcq, Jqq is nega- 
tive (ferromagnetic) and strongly dependent on Dcq as 
well as Dqq; this contrasts with the two-electron case 
[M] , where Jqq should be always positive (antiferromag- 
netic). This illustrates that even in the ground state the 
shallow-donor control electron can cause a significant per- 
turbation to the originally relatively weak exchange cou- 
plings in the two-qubit system. 




FIG. 4: (Colour online.) The probability 1 — ai of finding the 
control donor electron on the qubit cites for the case where a 
shallow-donor control and deep-donor qubits {i/ = 0.7) are in 
the spatial configuration shown in Figure [T] as a function of 
qubit donor distance Dqq and qubit-control distance Dcq- 
(a) control is in the ground state {v = 1.0); (b) control in the 
excited state {v = 2.0). Dqq varies from 1.0 to 12.0ao; H 
varies from to Wag. 



The transitions between different regimes can be 
largely understood from Figure [4] which shows the vari- 



ation of 1 



and hence of the probability (in the 



MuUiken sense) that the control electron is located on 
the qubit sites. We see that |ap ~ 1 throughout, except 
in the excited state when the control is relatively close 
to the qubits. In that case |ap <C 1 and the control 
state hybridizes strongly with the virtual orbitals on the 
qubit atoms. In this case the physics of the exchange is 
dominated by local interactions between the control and 
qubit spins similar to those in the ls2s excited state of 
He, where there is a ferromagnetic exchange splitting of 
6,422cm-i or 0.80 eV [37J. 

We focus on how the control electron alters Jqq and 
plot the ratio of this quantity with and without the con- 
trol present in Figure [5] As expected, when the control 
is very far away from qubits {Dqq 3> Dqq), the ratio is 
close to unity. When this condition is not satisfied (i.e., 
for Dcq ^ Dqq), Jqq is strongly modified and (over a 
large part of the region) ferromagnetic. In the excited 
state (Figure pTb)) the coupling is enhanced and ferro- 
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FIG. 5: (Colour online.) The ratio (logarithmic scale) of 
the qubit-qubit exchange couplings Jqq for two deep-donor 
qubits (i/ — 0.7) with and without a shallow-donor control. 
The spatial configuration is as shown in Figure IT] results are 
shown as a function of the qubit-qubit distance Dqq and 
the qubit-control distance Dcq- In (a) the control is in the 
ground state {f = 1.0); in (b) the control in the excited state 
{v = 2.0). The sign of the ratio is coded in color; red: pos- 
itive (antiferromagnetic interaction in presence of control), 
and green: negative (ferromagnetic). In the absence of the 
control, the interaction is always antiferromagnetic. 



FIG. 6: (Colour online.) The ratios J'' /J^ (logarithmic scale) 
of the exchange interactions for the control excited {J", with 
v — 2.0) to control-unexcited (J^, with v — 1.0) for a shallow- 
donor control and deep-donor qubits: (a) control-qubit ex- 
change JcQ, (b) qubit-qubit exchange Jqq. The sign of the 
ratio is coded in color: red: positive and green: negative. 



IV. RING-EXCHANGE MODEL 



A. Multi-center ring exchange 



magnetic almost everywhere, unless Dcq ^ Dqq. We 
will show in §IV| that the ferromagnetic interactions are 
a consequence of ring exchange among the three atoms. 

Finally, in Figure |6] we focus on the changes brought 
about by excitation (the 'on/off ratios' of the exchange). 
As would be expected, excitation almost always increases 
the magnitude of the exchange; Figure l6ja) (for Jqq) 
shows that this effect is largest for large Dqq and mod- 
erate Dcq, although there is a relatively complicated be- 
haviour of the sign of the ratio as a result of the differing 
crossovers between ferromagnetic and antiferromagnetic 
exchange in the ground and excited states (Figure |3|. 
FigureloFb) (for Jcq) shows a simpler pattern: the on/off 
ratio increases with Dcq and remains only weakly depen- 
dent on Don. 



In this section, we present an alternative calculation 
of the exchange constants by using Green's-function per- 
turbation theory. This approach has the advantage that 
it gives a clear picture of the origins of different contri- 
butions to the exchange. In particular it allows us to 
understand the origin of the ferromagnetism of Jqq in 
pil| owing to the presence of the shallow control atom, 
multi-body ring exchange process |19j become possible 
and higher-order ring exchange process (which can be fer- 
romagnetic) may dominate over the direct second-order 
exchange between the qubits provided Dcq <C Dqq. 

The essential physics of our system is similar to that 
in solid '^He [19 , despite the very different chemical na- 
ture of the material. In solid ''He, two-body exchange 
always leads to a Heisenberg Hamiltonian with anti- 
ferromagnetic exchange (Figure uh), while three-body 
ring exchange (Figure Up) and four-body ring exchange 
(Figure uh) lead respectively to ferromagnetic and anti- 



ferromagnetic contributions. 
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FIG. 7: Schematic view of possible exchange couphngs be- 
tween helium atoms in solid "^He: (a) conventional two- 
body exchange (antiferromagnetic); (b) three-body ring ex- 
change (ferromagnetic); (c) Four-body ring exchange (anti- 
ferromagnetic). 

As in the solid ^He system we need to include third- 
order and fourth-order processes (Figure [s]) in our 
Green's function perturbation theory calculation in or- 
der to capture these effects of control atom. 
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FIG. 8; Schematic view of multi-centre exchange processes in 
the control-qubit system: (a) third-order exchange scheme; 
(b) fourth-order exchange scheme. 



B. Green's functions time- independent 
perturbation theory 

Suppose the Hamiltonian for a system can be separated 
into an unperturbed part of Hq, and a perturbation V; 
H = Hq + V. We assume the eigenvalues and cigenfunc- 
tions of Hq are easily obtained. Green's functions Go{z) 
and G{z) corresponding to Hq and H, respectively, are 
Go{z) = {z- Ho)-\ and G{z) = (z - H)~\ G and Gq 
are connected by the Dyson equation 

G{z) = [\-G^{z)V]-^Ga{z) (17) 

= Go -b GoFGo + GoV^GqFGo -b . . . , (18) 

and the effective Hamiltonian which can mix a particular 
subset S of the eigenstates of Hq is 



A7Jeff(z) = PVP+PVQGo{z)QVP 

+PVQGo{z)QVQGo{z)QVP 



where Q = 1— P, P is the projection operator onto the set 
iS, and the 1st-, 2nd-, 3rd-order... energy shifts are ad- 
umbrated on the right-hand side of equation (fTOJ) . In this 



formula, AH^g depends on z, which should therefore be 
chosen to correspond to the energy of the states S. In the 
next section, we will show how we use Green's function 
perturbation theory to extract the exchange constants. 



C. Simple model of the three-centre problem 

From Figure |4] we see that the control electron resides 
on the control atom over the majority of parameter space; 
nevertheless. Figure [3] shows Jqq is ferromagnetic for 
much of this region. Therefore for simplicity we ignore 
the donor excited states, retaining only one orbital per 
site. We consider as parts of Hq all processes which do 
not transfer charge between centres (both single-electron 
terms and Coulomb interactions within one atom, and 
Coulomb interactions between different centres); the per- 
turbation V then contains hopping operators that trans- 
fer electrons between centres. Equation (17 1 then has a 



simple interpretation reflected in ( 19 1 : Go gives the prop- 



agation of uncoupled atoms, while V exchanges electrons 
between those atoms. 

Then we carry out our calculation in the following 
steps. First, we select a particular eigenvalue of the z- 
component of the total spin Sz (in the following discus- 
sion, we choose 5*2 = 2) ('^'z is a good quantum number 
since we ignore the spin-orbit coupling). Next we enu- 
merate all the possible configurations for this Sz- (For 
Sz = ^ the nine configurations are shown in Figure (|9|.) 
We choose which of these states we will retain as our 
set iS; typically, this is the low-energy subspace with one 
electron per donor, so that Coulomb intreactions are min- 
imized. Next we construct the matrices of Go and V 
based on all these possible states (see j jIV B I , and use 
these matrices to derive 2nd-, 3rd-, and 4th-order terms 
in Ai/eff according to equation ([T9|. (We suppose the 
higher terms are negligible.) Finally, as in the variatonal 
calculation (QIC I, we identify the values Jcq and Jqq 
by comparing the effective Hamiltonian in the last step 
with the Heisenberg Hamiltonian (15 1. 



.(19) 



D. Perturbation theory calculation 

1. The unperturbed hamiltonian Hq. 

We take for Hq a simple extended Hubbard-type 
Hamiltonian: 

Hq ^ Ecnc + EQ{nA + nB) 

+Ucnc^hci + UQifiATfiAi + fiB^fiBi) 
+VQQnAnB + VcQnc{nA + nB)- (20) 

Here Ec and Eq are the single-particle energies of the 
control and qubit atoms, Uc and Uq are on-site Coulomb 



|1) — C^-|-C(-,|Cgj^|U) \l) — C,^|C^^C^ j^|U) |o) — C^|C^|C|^| |U) 



|4) — c^|C(-,|Cg||0) |5) — Cp|C(-.|C^^ |0) |6) — c^|C^j^Cg-|-|0) 



|7) = c\^cl^cl^\Q) 18) = ct,^ct,^cl3jO) 19) = c^t4t4J0) 



FIG. 9: Schematic of nine possible states according to the 
calculation procedure in §IVC[ We define second-quantifying 



operators: 



CAo 



CCa 



cbct, where cr is j or |, 



and c\^ is electron-creating operator on the state of qubit A, 
etc. 

interactions in the control and qubit, and Vqq and Vcq 
are respectively Coulomb interactions between electrons 
in the two qubit sites, and between electrons in the con- 
trol and qubit sites, fiia- is the number operator for elec- 
trons with spin a on site i, and hi = hi^ -t- fiii. 

We take for the perturbation V the hopping terms that 
transfer electrons respectively between the control atom 
and the qubits, and from qubit to qubit. If we assume 
the corresponding amplitudes tcQ , iqq are real, then we 
can write 



V 



st 



st 






[tcqi&Aa^Ca + &Ca^Aa + c\^CCa + &Ca^Ba) 



+tQQ{c\„CBa + C^Ba^Aa)]. (21) 

where Cia annihilates an electron with spin <j at site i. 



2. Perturbation theory and reproduction of the Heisenherg 
spin Hamiltonian 

Now we have Go and V, we can reproduce the Heisen- 
berg Hamiltonian from the 2nd-, 3rd- and 4th-order per- 
turbations to find the corresponding contributions to 
Jqq and JcQ according the calculation procedure in 



Jqq 
JcQ - 



J, 



(2) 



J 



QQ 

(2) 
CQ 



7(3) 

'-^QQ 

'^CQ ^ 



.(4) 

-^QQ 

t(4) 
'^CQ- 



(22) 
(23) 



r(») t(») 



■^QQ' ^CQ^ (* = 2,3,4) are functions of Vqq, Vcq, Uq, 
Uc, tQQi and tcq', all the conributions are given explic- 



itly in Appendix [Bj 

Motivated by the known physics of the defect problem, 
we choose suitable values for these variables to illustrate 
our results. We expect the hopping integrals Iqq and tcQ 
to decay exponentially with Dqq and Dcq respectively, 
at rates determined by the sum of the orbital exponents 
of the relevant atom pairs. We also expect Vqq, Vcq < 
Uc < Uq because the on-site Coulomb interaction should 
be larger than the Coulomb interaction between different 
centres, and the qubit spatial state is more localized than 
the control state. 

For definiteness we take the following values for the 
parameters: Uq = 1.0, Uq = 1.5, Vqq = 1/{2Dqq), 



Vcq - 1/{2Dcq), 



QQ - 
= 0.7, 



-10- 
■r. - 1.0. 



Using these pa- 



rameters, we calculate Jcq^Jqq shown in Figure (10 1 




FIG. 10: (Colour online.) The exchange couplings (a) Jcq 
and (b) Jqq calculated by perturbation theory as a function 
of qubit distance Dqq and qubit-control distance Dqq- The 
sign of exchange coupling is coded in color; red: positive (an- 
tiferromagnetic) and green: negative (ferromagnetic). Dqq 
varies from 1.0 to 12.0; H varies from to 10. 

Figure [lO] shows us that when the qubit-qubit trans- 
fer amplitude Iqq is much weaker than the control-qubit 
transfer tcQ, ring-exchange will dominate over simple 
two-centre exchange and hence lead to ferromagnetic ex- 
change couplings. The variations of exchange coupling 
shown in Figure 10 are qualitatively similar to those 
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shown in Figure l2{a) and Figure l3|a): when Dcq > 
Dqq, the exchange coupling is antuerromagnetic, while 
for some configurations where Dcq « Dqq the ex- 
change coupling is ferromagnetic. Note that in the values 
assumed for the parameters, the prefactor for Iqq is cho- 
sen so that it is always much smaller than tcQ] we find 
that the magnitude of this prefactor tunes the bound- 
ary between the ferromagnetic and antiferromagnetic re- 
gions, as might be expected because of its effect on the 
relative magnitudes of third- and fourth-order processes. 
For the excited control, tcQ is larger and the ferromag- 
netic region is consequently further extended. Ultimately 
with these parameters a negative (ferromagnetic) part in 
the qubit-qubit exchange coupling Jqq arises at interme- 
diate qubit-control separations. 



V. STATISTICAL DISTRIBUTION OF 
EXCHANGE INTERACTIONS 

Since there is a strong configuration dependence of 
the exchange interactions it is natural to ask about 
their probability distributions at given impurity con- 
centrations, since these distributions will determine the 
system's response to macroscopic, spatially averaging, 
probes. We proceed as follows: we define a cluster 
of three spins to consist of a qubit atom, its nearest- 
neighbour qubit atom, and the control atom nearest to 
their mid-point. Then, assuming the qubits and controls 
are independently and uniformly distributed (i.e. ne- 
glecting both defect-defect correlations and the unerlying 
lattice structure), the probability P{Dqq, H) cIDqq dH 
that the qubit-qubit distance is between Dqq and Dqq + 
cLDqq while the control's distance from the mid-point is 
between H and H + dH, is given by a simple generaliza- 
tion of the argument of Chandrasekhar f35j as 

P{Dqq,H) = 16vr2i?^Qif2^Qnce[-*("«^?^«+"-^^)], 

. (24) 
where uq and nc are the volume densities of qubit and 
control atoms respectively. For purely exponential in- 
teractions between defect pairs, the averages over this 
type of distribution may be performed analytically HO] ; 
since we know that in this case the interactions depend 
on both Dqq and H, and are not purely exponential, we 
must resort to a numerical average. 

Our key approximation is to assume that the exchange 
couplings computed in ^11 (i.e., in an isosceles configu- 
ration with both qubit-control distances equal) are rep- 
resentative of all configurations with the same values of 
Dqq and H; we suspect that this approximation will 
lead to excessive delocalization of the control electron 
and hence to an over-estimate of the exchange. We com- 
pute the induced distribution of exchange couplings from 
our previous results by summing over the same mesh 
used to generate Figures [2] and [3] accumulating a suit- 
ably weighted histogram of exchange values by convolut- 
ing the delta-function contributions from each calculated 



point with a Guassian whose width is chosen to produce a 
smooth curve without introducing excessive broadening. 
At the densities chosen, this mesh accounts for approx- 
imately 98% of the normalization of P{Dqq,H). The 
results are shown in Figures (111 and (IT2| , in the form of 



distribution functions for log \J\ with the contributions 
from positive and negative J (i.e. antiferromagnetic and 
feroommagnetic parts) displayed separately. From Fig- 
ure fll|a), we can see that when the control is in the 
ground state the distribution of qubit-control exchange is 
overwhelmingly antiferromagnetic and qualitatively sim- 
ilar to those calculated for defect pairs of similar binding 
energies at comparable densities [17! . However, when the 
control is excited (Figure ITlJb)) the antiferromagnetic 
distribution narrows dramatically and shifts to higher 
I J|, both because of the slower decay with distance and 
because of the significant ferromagnetic region discussed 
in §III[ At the same time the ferromagnetic part to the 
distribution grows. 

Similarly, when the control is in the ground state the 
distribution of Jqq (Figure 12) is predominantly antifer- 
romagnetic and very close to that for deep-donor pairs 
at the same density [17j. The small ferromagnetic part 
is due to the ring exchange described in |IV[ When the 
control is excited, this ferromagnetic region is greatly en- 
larged owing to the enhanced tcQ , which leads to ring ex- 
change dominating the two-donor exchange as illustrated 
in 3171 



VI. DISCUSSION AND CONCLUSION 

We have described a variational method for calculat- 
ing the exchange interactions between three electrons in a 
control-qubit system of the type envisaged for optically- 
controlled quantum information processing in silicon [S\. 
We have deliberately kept the model simple, yet rich 
enough to capture the coupling of the changing electron 
distribution to the spin configuration as the geometry 
and degree of excitation are altered. 

The principal conclusion is that a three-centre system 
such as this is significantly different from the sum of two- 
electron exchange interactions that would arise from con- 
sidering each defect pair separately. In particular the 
excitation of the control electron can alter the 'qubit- 
qubit' coupling Jqq as well as the 'control-qubit' cou- 
pling JcQ- This is partly because of the delocahzation 
of the electronic states across all three centres (partic- 
ularly when the control electron is in an excited state), 
and partly because of the presence of three-centre ring- 
exchange processes (which we have also described per- 
turbatively) . Most strikingly, we predict that either or 
both the control-qubit coupling Jqq and the qubit-qubit 
coupling Jqq can become ferromagnetic; this can happen 
through either of the above mechanisms. The excitation 
of the control greatly enhances both the ferromagnetism 
and the degree of charge transfer; however, it is interest- 
ing that Jqq can be ferromagnetic even when the control 
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FIG. 11: (Colour online.) The probability distribution of 
log I JcQ I for the control-qubit coupling, shown as a function 
of JcQ (on a logarithmic scale) for the control atom in (a) its 
ground state and (b) its excited state. Red solid curve: dis- 
tribution of antiferromagnetic couplings; green dashed curve: 
ferromagnetic couplings. The density of qubit atoms was 
taken as 0.00396 (aj)"^ = 5.42 x 10"cm"^, and the den- 
sity of control atoms as 0.00198 (aS)~^ = 2.71 x 10"cm"^. 
The plot was constructed by sampling over the 111 x 101- 
point grid used to create Figure [2J smoothed by convolving 
each point with a Gaussian in log \Jcq\ of standard deviation 
0.224. 



is in the ground state. The quahtative agreement be- 
tween our variational results and those of the simple per- 
turbation theory (which includes only the ring-exchange 
mechanism) shows that the latter is especially important 
when the control is in the ground state. The ferromag- 
netic couplings we predict could be exploited for the ini- 
tialization or manipulation of donor spins. 

As might be expected, when the control electron moves 
further away, Jqq returns to the value expected for a 
deep donor pair without any control species |17j . This is 
also reflected in the statistical distributions of exchange 
values, where the distribution of Jqq with the control 
in the ground state is quite similar to the exchange dis- 
tribution for deep-donor pairs calculated in [17]. How- 
ever, when the control electron is located between the 
two qubits the magnitude of Jqq increases significantly, 
owing to the mediation by the control electron; this is es- 



(b) 

FIG. 12: (Colour online.) The probability distribution of 
log I Jqq I for the qubit-qubit coupling as a function of Jqq 
(on logarithmic scale) for the control atom in (a) the ground 
state and (b) the excited state. Red solid curve: distribution 
of antiferromagnetic couplings; green dashed curve: ferromag- 
netic coupling distribution. The parameters used to construct 
the plot were as for Figure [TTj 



pecially significant when the qubit-qubit distance is large. 
In this limit Jqq can also become ferromagnetic, because 
the fourth-order terms containing tcQ which give a ferro- 
magnetic contribution dominate. When the qubit-qubit 
distance is relatively small, the second-order terms in- 
volving Iqq dominate Jqq, and the control electron is a 
minor perturbation on the two-qubit sub-system. 

Our results show that the process of excitation and 
de-excitation can produce substantial but quite complex 
changes in the exchange couplings and charge distribu- 
tions in these model defect complexes. They provide a 
theoretical demonstration of the switching of exchange 
couplings by excitation; a quantitative understanding of 
this process will have implications both for optically- 
controlled magnetism in silicon, and for quantum com- 
putation. 
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APPENDIX A: COMPUTING THE 
COEFFICIENTS IN THE SPIN HAMILTONIAN 



Equation (15) is the most general rotationally- 



invariant three-spin Hamiltonian that is even under time 
reversal. Total spin angular momentum is a good quan- 
tum number, and standard angular momentum coupling 
schemes |38| show that three spin- 1/2 particles can be 
coupled to a quartet \Q) — {[5, 5] 175)3, and to two 

different doublets, \Dl) = {[s'sIl'sIi distinguished 



by the intermediate angular momentum L, which can be 
or 1. 



Since the full spin Hamiltonian (15 1 is traceless, we 



may conveniently write its eigenvalues as 



5d„ = 


-e - (5/2 (doublet) 


5cl = 


-e -1- 6/2 (doublet) 


Eq = 


e (quartet) 



(Al) 



It is sufficient to work with the three states having to- 
tal spin component 1/2 (|001), |010) and |100)); they are 
only coupled amongst themselves, and within this sub- 
space the Hamiltonian is 



JcQB - JcQA 
"^Jqq 

"^JCQA 



which has one eigenvalue 



2^QQ 2JcQA 

'JcQB + JcQA — JqQ "^JcQB 

"^JCQB —JcQB — JcQA + JqQ 



JqQ + JcQA + JcQB 



= Ec 



(A2) 



(A3) 



(corresponding to the component of the quartet with magnetic quantum number=l/2), while the other two eigenvalues 
are 



1 



e ± 7.^ {J COB + JcQA + Jqq) - JcqaJqq - JqqJcqb - JcqbJcqa 



(A4) 



(corresponding to the two doublet states). 



Now, in the symmetric case Jcqa = Jcqb = JcQ, it 
is straightforward to show that 

C^{{sa + sbY) (A5) 

commutes with the Hamiltonian. Since C is the operator 
which corresponds to the intermediate angular momen- 
tum L described above, we may characterize the two dou- 
blet states by their intermediate L value. Once more it is 
convenient to work within the spin 1/2 subspace, where a 
short calculation shows that the doublet state with L = 
has the eigenvector {1, — 1,0}-^/\/2, with eigenvalue 

^, (A6) 

1 has eigen- 



Edo - 



whereas the doublet characterized by L 
value 



Edi — 



J, 



CQ 



(A7) 



with eigenvector {1, 1, — 2}^/\/6. The quartet state, of 
course has eigenvalue 

JcQ . Jqq 



Er 



(A8) 



and eigenvector {1, 1, Ij-^/'s/S. Thus, the simplest way 
of determining the exchange parameters Jqq and Jcq is 
to identify the doublet state with L = and the quartet 
state and use the relations (16 1 given in the main text. 



Notice it is important to identify the doublet states 
by their wavefunction; their energies alone are not suffi- 
cient. The eigenvalues of the Hamiltonian defined by the 
exchange parameters Jcq and Jqq are identical to those 
of the Hamiltonian obtained by the substitutions 



Jcq ^ o("^c'Q + 2Jqq); 
Jqq — > ^{"^JcQ - Jqq)- 



(A9) 



However, the eigenweciors of the doublets interchange 
their properties under this substitution. This is related 
to the sign ambiguity in the square root in equation A3, 
which makes the sign of (5 in A 1 ambiguous. 



12 



APPENDIX B: ANALYTIC RESULTS FROM 
PERTURBATION THEORY 



The lowest-order (second-order) contributions to the exchange are 



7(2) _ 



-Uq + 2VcQ - Vqq 



9/2 



Uc'' 



(Bl) 



7(2) 



iti 



Vqq - Uq 



(B2) 



The third-order contributions are 



J, 



(3) _ 



QQ 



while the fourth-order terms are 



^(3) ^tcQ^QQ 


'^^CQ^QQ 




^'^ {VQQ-Ucy { 


-Uq + 2Vcq-VqqY' 




'^^QQ^CQ 


^^QQ^CQ 




{-Uq + 2Vcq-Vqq){Vqq- 


- Uq) {Vqq - Uc){Vqq - 


-Uq) 


'^^QQ^CQ 


'^^QQ^CQ 





(B3) 



(B4) 



{-Uq + 2Vcq - Vqq)2 {Vqq - Uc)^ 



'cq 



^tcQ 



^'^CQ 



{-Uq + 2VcQ - Vqq){Vqq - Uc){Vqq - Uq) {-Uq + 2VcQ - Vqq)^{Vqq - Uq) 



(B5) 



"^^CQ 



r,.2 f2 

'^^QQ'-CQ 



r,.2 f2 

'^''QQ^CQ 



{Vqq - Uc)^{Vqq - Uq) {-Uq + 2VcQ - Vqq)^ {Vqq - Uc)' ' 






^t%Q 



^^QQ^CQ 



{-Uq + 2VcQ - Vqq){Vqq - Uc){Vqq - Uq) {-Uq + 2VcQ - Vqq)^Vqq - Uq) 



(B6) 



4+2 f2 

^'-qq'^cq 



^^QQ^CQ 



4/2 +2 

^'^qq'-cq 



{Vqq - Uc)^{Vqq - Uq) {-Uq + 2Vcq - Vqq){Vqq - Uq)^ {Vqq - Uc){Vqq - UqY' 
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